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THE  EXISTENCE  AND  CALCULATION  OF  SOLUTIONS  OF  CERTAIN 
INTEGRO-DIFFERENTIAL  EQUATIONS  IN  SEVERAL  DIMENSIONS 

INTRODUCTION 

In  this  paper,  boundary  and  initial  value  problems  will  be  solved  for 
linear  integro-differential  equations,  in  several  dimensions,  of  a  type  in¬ 
cluding  linearized  forms  of  Boltzmann's  equation.  The  methods  of  solution 
are  based  on  finite  differences  and  are  thus  suitable  for  calculation. 

These  methods  are  indirect,  for  our  difference  scheme  is  proved  to 
converge,  not  for  the  original  problem,  but  for  a  modification  of  it,  whose 
solution  approximates  that  of  the  original.  A  better  difference  scheme,  or 
a  better  proof  Of  convergence,  might  eliminate  the  need  to  modify  the  original 
problem.  This  need  to  modify  the  problem,  however,  may  be  real,  desirable 
in  numerical  practice  as  well  as  in  theory. 

Most  discussions  of  the  integro-differential  equations  treated  here 
occur  in  the  context  of  neutron  transport  theory,  of  which  a  comprehensive 
account  is  to  be  had  in  B.  Davison's  book  [l].  Calculation  techniques  in 
transport  theory  are  also  surveyed  in  R.  D.  Richtmyer's  book  [2]  and  in  a 
subsequent  paper  [3]  of  E.  H.  Bareiss.  The  first  proof  that  general  multi¬ 
dimensional  problems  resembling  ours  have  solutions  was  given  by  semi— group 
methods  by  K.  JBrgens  [a].  We  do  not  refer  here  to  the  extensive  literature 
on  problems  in  one  or  two  dimensions. 
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1.  Statement  of  problem.  Notation.  Main  Results.  Let 
2&  :  lxrl  <  c>  r  =  1,  ...»  d, 

designate  a  d-dimensional  cube  of  edge  2C,  the  symbol  x=(x^,Xg,  ...,x^) 
referring  to  any  of  its  points.  Let  y  =  (y^>  •••»3r(j)  an<^  y'  = 
denote  any  d-dimensional  vectors  and  dy'  =  dyj^  . ..dy^  the  d-dimensional 
element  of  volume,  and  let  t  be  an  independent  variable  on  the  semi— axis 
t  >  0.  The  equations  considered  are  those  of  the  form 

1.1  u .  '+  2  y  u  +  c(x,y,t)u  =  \K(x,y,t,y')  u(x,y' ,t)dy'  +  g(x,y,t) 

1  r=l  r  xr  J 

with  coefficient  c(x,y,t)  and  inhomogeneous  part  g(x,y,t)  defined  in  a 
(2n+l)— dimensional  cylindrical  drum, 

ST:  x  tjJ,  0  <  t  <  T,  -oo  <  yg  <  oo  (T=const.  >  0),  s  =  l,...,d 

and  with  kernel  K(x,y,t,y' )  defined  in  a  (3n+l)-dimensional  cylinder 
over  ST, 

T..!  (x,y,t)e  ST,  -oo<  y*  <  oo,  s  =  1,  ...,  d. 

T  is 

The  domain  of  integration  here  has  been  taken  as  the  entire  d-dimensional 
y '-space,  although  a  bounded  portion,  or  a  lower  dimensional  subset,  would 
have  been  equally  acceptable.  Unlike  Jorgens,  we  do  not  exclude  a  neighbor¬ 
hood  of  the  origin  from  the  domain  of  integration. 

Later  on,  S  will  denote  the  base,  t  =  0,  of  S  *.  will  be 

o  1 

the  set  of  points  of  S^,  on  any  of  the  planes  yr  —  0,  r  =  1, ...,  d.  For 
e  >  0,  will  be  an  e-neighborhood  of  in  ST. 
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4 


will  denote  an  (n+l)-dimensional  cylindrical  drum  with  base  £  , 

namely 

XT:  x  e  %  ,  0  <  t  <  T  . 

By  %  will  be  meant  the  set  of  points  of  ^  on  any'  of  the  planes 

yr  -  0,  r  =  1,  ...,  d  . 

The  initial  condition  imposed  is  of  the  form 
1.2  u(x,y,0)  =  <j>(x,y)  . 

The  boundary  condition  is 


1.3 


u(x,y,t)  =  0  when  z  y.N  (x)  <  0  , 

J  J 


where  (N.(x))  denotes  the  outward  normal  to  £  at  x,  a  point  of  the  boundary 
J 

not  on  an  edge.  This  boundary  condition  can  be  written  as 

1.3'  u(x,y,t)  *  0  for  x  e  B  , 

r»y 

where 


B 


r,y 


|xsl  <  C,  s  4  r 


x  =  C  =  C,  if  y  <  0 
r  r,y  ’  •'r 


=  — C,  if  y  >0 


is  an  appropriate  one  of  the  two  faces 

Br:  xr  =  C>  lxsl  <  C  for  3  4  r  , 


B  :  x  =-C, 
r  r  ’ 


<  C  for  s  4  r  , 


of  %  normal  to  the  x^-axis .  When  yr  =  0,  no  boundary  condition  is  imposed, 
and  Br  in  this  case  is  to  be  taken  as  empty.  For  ^ 

B  (r=l, ...,d),  d  mutually  perpendicular  planes  in  d-dimensional  space, 

r»y 


3 
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all  Intersect  in  one  point. 

By 


B1  : 
r,y 


/'IxJ  <  C,  xs  =f  Cs>y  for  s  4  r 


x  =  C  if  y  >  0 
r  jt 

=  -C  if  yr  <  0 
=  ±C  if  yr  =  0 


we  shall  denote  the  face  of  Jp  normal  to  the  xr-axis  on  which  u  goes  un¬ 
prescribed. 

U  will  denote  the  union  of  the  planes  y  =  0,  s  =  1, . . . ,  d;  the 

s 

equation  of  U  is  y^**^  =  0- 

We  shall  seek,  depending  on  the  precise  hypotheses,  either  "weak"  or 
"strong"  solutions  of  1.1,  1.2,  1.3. 

DEFINITION  1:  By  a  "strong"  solution  in  ST  will  be  meant  a  bounded 
function  u(x,y,t)  that,  for  constant  x,t,  is  measurable  with  respect  to 
y,  for  constant  y  ^  U  is  absolutely  continuous  with  respect  to  t,x  in 
T  -  fulfills  the  initial  and  boundary  conditions  1.2  and  1.3, 

and  satisfies  1.1  at  almost  all  points  of  S^. 

"Weak"  solutions  are  functions  that  satisfy  the  required  conditions  in 
a  certain  integral  sense.  To  arrive  at  an  appropriate  formulation,  imagine 
u  to  satisfy  1.1,  1.2,  and  1.3  strictly,  and  this  and  the  other  functions 
involved  to  admit  the  various  operations  we  shall  perform.  We  first  multiply 
equation  1.1  by  <p(x,y,-t)F’  (u),  where  F(u)  and  <p(x,y,t)  are  continu¬ 
ously  differentiable  functions  for  all  values  of  their  arguments.  Then  we 
integrate  over  an  arbitrary  polyhedral  domain  D  contained  in  the  drum 
Letting  D  denote  the  boundary  of  D  and  (v.,v.)  ,  the  outward 

2.  T i  1—  ±9  •  •  •  f  0. 
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normal  to  this  boundary,  equation  1.1  becomes,  upon  integration  by  parts, 
^  <P  F(u)(vt  +  2  yr  vr)dS  -  J' F(u)(<pt  +  2yr  <px  )dxdt 


-J  <p  F*(u)  -cu  +  ^Kudy1  +  g| 


dxdt  . 


Now  we  specialize  D  to  0  <  x  <  T,  take  F(u)  =  u,  and  impose  the 

condition 

1.5  ®  =  0  on  B'  ,  r  =  1,  ...,  d. 

In  view  of  1.3*  and  1.2,  the  integral  over  D  in  1.4  thereby  becomes 

I®  udx|t-T  “  J  ®  ^  * 

and  relation  1.4 

=  JV  ■t’dx  +  J*u(<pt  +  2yp<Px  )dxdt 

,,  r  X  t 


+  j  <p(— cu  +  J* Rudy'  +  g)dxdt  . 


We  thus  arrive  at 


DEFINITION  2:  u(x,y,t)  is  called  a  "weak"  solution  of  1.1,  1.2,  1.3 
if  relation  1.6  holds  for  all  vectors  y,  all  x  in  the  interval  (0,T), 
and  all  functions  <p(x,y,t)  continuously  differentiable  in  ST  that  satisfy 


Appropriate  function  classes  for  u  and  the  other  quantities  involved 
will  be  indicated  below. 
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Under  rather  light  hypotheses,  e.g.,  hypotheses  (i)-(iv)  below,  a 
bounded,  absolutely  continuous  (in  )  weak  solution  is  a  strong 

solution.  This  is  so  because  the  procedure  that  had  led  formally  to  1.6 
can,  in  the  case  of  an  absolutely  continuous  u,  be  reversed.  Conversely, 
a  strong  solution  is  weak. 

Let  V(ST)  denote  the  space  of  bounded  functions  v(x,y,t)  Lipschita- 
continuous  on  S,p.  Define  V'(S^)  as  its  completion  under  the  norm 


M*  -  SUP 


OitiT 
all  y 


|  J(v(x,y,t))2dxj 


3/2 


define  V" (S^,)  as  the  completion  of  V(S^)  under  the  second  norm 


|v  [  I "  =  sup1 


ijLj*  (v(x,y,t))2dxdt|  j 


let  V(Sq)  denote  the  space  of  bounded  functions  (j)(x,y)  Lipschita- 

continuous  on  S  and  V'(S  )  the  completion  of  this  space  under  the  norm 
00 

M<j>||'  =  sup  \  \  (<J>(x,y))2dxl 


y  \  *  ) 

Any  member  of  V *  ( S^, )  is  a  function  w(x,y,t)  such  that 

3/2 


w(y. 


r,t)  =  |  J(w(x,y,t))2dx^ 
is  continuous  in  y,t;  similar  statements  hold  respecting  the  members  of 


Vn(S^.)  and  V'(So)  .  To  see  this,  let  v^(x,y,t),  k=  1,2,...  , 
belonging  to  V(ST),  approximate  w  in  V'(ST)  .  Thus, 


lim 


sup 

k,m-»  00  OltlT  J  * 
ally  t 


1' 


vm(x,y,t))  dx  =  0  . 
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This  implies  that  the  continuous  functions 

1/2 


Wk(y,t)  =  |  J(vk(x,y,t))2dx| 

converge  uniformly  with  respect  to  y,t  ,  since  by  the  triangle  inequality 


we  have 


|Wk(y,t)  -  Wm(y,t)|  <  |^J(vk(x,y,t)  -  va(x,y,t)  )2dx| 


1/2 


and,  therefore, 


lim  sup  |W,  (y,t)  -  W  (y,t)|  =  0  . 
k,m-*  00  0it4T  K  m 

all  y 

Hence.  W  =  lim  W,  is  continuous,  as  asserted. 

’  k 

We  shall  prove  the  existence  of  weak  solutions  of  our  problem  under  the 
following  hypotheses: 

(i)  The  coefficient  c(x,y,t)  is  bounded  and  belongs  to  V"(ST). 

(ii)  The  inhomogeneous  part  g(x,y,t)  belongs  to  V"(ST). 

(iii) The  kernel  K(x,y,t,y' )  satisfies  the  following  conditions: 

(a)  K  >  0  . 

(b)  A  sequence  of  non-negative , continuously  differentiable  f  amo¬ 
tions  Kjxjyjtjy' )  ,  m=  1,2,...,  exists  such  that 

% 

|K(x,y,t,y' )  -  Km(x,y,t,y* ) |dy'  <  em 


!< 


in  S^,  ,  where  the  are  constants  that  approaGh  zero  as  m«^oo. 

(c)  K(x,y,t,y« )  <  Ko(y,y' )  ,  where  KQ(y,y' )  is  integrable 

with  respect  to  y'  and  also  integrable  with  respect  to  y  on  any  sphere 
|y|  <  const.,  and  where 
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jKo(y,y')dy'  <  kQ  , 


kQ  being  a  constant. 


(d)  For  all  y  , 


iim  J^J|K(x,y,t,y')  -  K(x,y1,t,y' )  |dyj^  dxdt  =  0 


(iv)  The  initial  data  function  (jl(x,y)  belongs  to  V'(So)  . 

If  u  is  in  V'(S^)  and  the  above  hypotheses  are  satisfied,  the 
integrals  in  1.6  will  be  continuous  with  respect  to  y,  t.  We  see  this 
first  for  the  integral  on  the  left,  which  can  be  expressed  as  a  difference, 

J<pudx  =  ^  J\<p  +  u)2dx  -  ^  J'Ctp  -  u)2dx  , 

of  two  terms  already  known  to  be  continuous.  Analogous  considerations  apply 
to  the  first  and  second  integrals  on  the  right  and  also  to  ^ cpcudxdt  and 
J*  <pgdxdt  .  The  integral 


Kvdy'  , 


in  which  v  =  cpu,  is  easily  seen  from  Schwarz'  inequality  and  hypothesis 
(iiic)  to  be  continuous  in  t  uniformly  with  respect  to  x,y  .  Hence,  all 
that  remains  is  to  prove  that 


(K(x,y,t,y')  -  K(x,y1,t,y'))v(x,y,,t)dy'^  dxdt 


tends  to  zero  as  y^— ►  y.  Fubini's  theorem  and  Schwarz'  "inequality  show  the 
latter  integral,  in  absolute  value,  to  be 
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i 

s 


l 


j 

i 

1 

: 


f 


f 

{ 


<  Jdy'  |^J|K(x,y,t,y'  )-K(x,y1,t,y' )  |dxdtj  j  ^JlKU.y^y' )  -  K^y^t,?' ) 


<  j  J|K(x,y,t,y' )  -  K(x,y1,t,y' )  |dy'^ 


(v(x,yf ,t))  dxdt 

} 


? 


2/2 


•)  2/2 

dxdt  V 


■ili"  x,y,t,y' 


)  -  K(x,y^,t,y' ) I  (v(x,y‘ ,t))2dxdtdy' 


2/2 


and  by  (iiic)  to  be 


2/2 

<  (2ko)  | | v | | 


"  j  J  ^J|K(x,y,t,y*)  -  K(x,y1,t,y' )  |dy^  dxdtj, 


This  tends  to  zero  by  (iiid).  Thus,  all  the  integrals  that  enter  1.6  are 
indeed  continuous  in  y,r,  as  asserted. 


THEOREM  1  (Existence  of  weak  solutions):  Under  hypotheses  (i)-(iv), 
the  problem  of  satisfying  1.1,  1.2,  1.3  has  a  weak  solution  belonging  to 

v'(sT)  . 

The  weak  solutions  we  shall  actually  construct  belong  to'  a  narrower 
class  than  that  ascribed  to  them  above,  and  in  the  narrower  class  they  are 
unique,  i.e.,  are  uniquely  determined  by  their  initial  data.  These  weak 
solutions,  namely,  all  have  the  property,  hitherto  unmentioned,  of  satisfy¬ 
ing  the  inequality 


1.7  J(u(x,y,x))2dx  <  f(<j>(x,y))2dx  +  2  f|-cu2  +  uJ*Kudy'  +  guj  dxdt 

X  ~£  x 

for  all  vectors  y  and  all  x  in  the  interval  (0,T).  This  inequality. 
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like  1.6,  is  formally  deducible  from  1.4  in  which  we  would  take 
D  =  *£^,F(u)  =  u2  and  f  =  1  to  obtain 

Ju2(vt  +2f  yrvy)dS  =  2  J'  |-eu2  +  u  f  Kudy’  +  gu*l  dxdt  . 

*■<  *t  5 

The  left  member  of  the  latter  relation,  however,  equals 
JuU.y.tl^dx  -  Ju(x,y,0))2dx  +  Zg  C  dt  J  u2  'J^T  d^x  , 

£  £  *  Bs 

d  x  denoting  the  element  of  area  on  B'  ,  while  on  B'  2  y  v  >  0.  From 
s  s  s  r  r  ■— 

these  remarks,  inequality  1.7  easily  follows,  at  least  for  strict  solutions 
of  suitably  regular  problems.  For  weak  solutions  as  described,  the  inequality 
will  be  seen  to  carry  over  by  closure. 

THEOREM  la:  The  weak  solutions  rendered  in  Theorem  1  constitute  a 
linear  class  W  each  of  whose  members  u(x,y,t)  is  subject  to  inequality 
1.7  with  appropriate  <|>(x,y)  and  g(x,.y,t)  . 

THEOREM  2  (Uniqueness  of  weak  solutions):  The  members  of  W  are 
uniquely  determined  by  their  associated  (J)'s  and  g's:  the  member  correspond- 
ing  to  (j)  =  0  and  g  =  0,  in  particular,  is  u  =  0. 

When  c,g,(}>,  and  K  are  requisitely  smooth  (continuous  differentiability 
would  be  enough),  and,  in  addition,  K  is  sufficiently  attenuated  near 
y'  =  eq  a  strong  solution  of  the  problem  exists.  Respecting  c,  g,  and  <J>, 
the  following  assumptions,  added  to  the  previous,  suffice: 

(i)Q  c(x,y,t)  and  g(x,y,t)  axe  uniformly  H&lder-continuous  with 
respect  to  x  ,  s  =  1, ...,d,  with  Hblder  exponent  a(0  <  a  <  l)  and  HtJlder 
constants  denoted  by  cq  and  bQ,  respectively. 
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(i)  c(x,y,t)  and  g(x,y,t)  are  uniformly  Lipschita-continuous  with 
respect  to  t. 

(iv) ^  4>(x*y )  is  Lipschita-continuous  with  respect  to  x  with  a 
uniform  constant  denoted  by  ^  . 

(it  might  be  expected  that  these  conditions  should  be  of  the  same  kind, 
all  Lipschitz,  or  all  Httlder,  conditions.  With  all  Lipschitz  conditions,  the 
solution,  however,  still  could  not  be  proved  to  be  more  than  Holder  con¬ 
tinuous  (with  any  exponent  <1  in  an  appropriate  region),  and  with  only 
Hfllder  conditions  nothing  could  be  proved  at  all.  Therefore,  the  three 
hypotheses  have  been  left  of  different  kinds.) 

Respecting  K,  four  new  assumptions  are  made: 

(v) ^  For  t  >  0,  At  >  0, 

I  dy*  <  kx  , 

where  k^  is  a  constant. 

(r) 

^2  ^*r  4  0»  x  =  (x^, . .  •,xjv_^,xr+  ^r,xr+^, . .  .fX^)  , 

we  have 

f  |K(x(r),y,t,y')  -  KU,?,t,?')|  dy,  <  k  f  r=l,...(d, 

V  I  Ax  [  a 

1  r  1 

where  k  is  a  constant, 

a 

(v)3  Let  yjrj  =  (y^,...,  y^,  y^, . .  .,y^) .  A  positive  X  and,  for 
each  r  =  1, ...,d,  a  function  Kr(y[r))>  exist  such  that,  if  Jy^ |  <  X.  , 

K(x,y,t,y')  <  Kr(yjrj) 

and 
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^Ir)  <  k 


for  r  =  1, where  k  is  a  constant. 

(v).  A  constant  K..  exists  such  that 
4 

J  K(x,y,t,y' )  |y^|  dy'  <  ^  ,  r=  l,...,d. 

We  can  now  state 

THEOREM  3  (Existence  and  uniqueness  of  strong  solutions):  Under 
hypotheses  (i)-(iv),  (i)^,  (i)t»  (iv)^,  (v)^  -  (v)^,  the  problem  of  satis¬ 
fying  1.1,  1.2,  1.3  has  a  strong  solution,  which  is  Holder-continuous  in 
S^,  -  S(o)  and,  for  any  e  >  0,  uniformly  so  in  S,p  -  S^  .  This  solution 


is  unique. 

The  paper  is  organized  as  follows: 

Section  1:  Statement  of  problem.  Notation.  Main  results. 

Weak  and  strong  solutions  are  defined  and  the  main  results  as  to  their 

existence  and  uniqueness  stated. 

Section  2:  Uniqueness  and  continuous  dependence.  Weak  solutions  of  the 
class  W  defined  above  (W  includes  all  strong  solutions)  are  proved 
to  depend  continuously  on  their  initial  data  and  thus,  in  particular, 
to  be  uniquely  determined  by  their  data  (Theorem  2) .  Since  the  weak 
solutions  obtained  (Theorems  1  and  la)  belong  to  W  ,  they  are  thus 
seen  to  be  independent  of  accidental  features,  like  selection  procedures, 
pertaining  to  the  construction  process. 

Section  3:  Positivity.  Monotonic  dependence  of  solution  on  data,  co¬ 
efficient,  kernel,  and  inhomogeneous  part  of  equation.  The  solution  is 
shown  to  be  positive  when  the  quantities  indicated  have  appropriate 
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fixed  signs.  All  solutions,  therefore,  depend  monotonically  upon  these 
quantities. 

Section  4:  Truncated  problems.  Proof  of  Theorems  1  and  la,  A  special  kind 
of  problem  is  defined  that  later  proves  solvable  by  finite  differences. 
Anticipating  that  the  solutions  of  these  special  problems  exist,  we 
here  show  how  they  can  be  used  to  approximate  the  solutions  of  arbitrary 
problems.  Thus,  Theorems  1  and  la  are  proved. 

Section  5:  Difference  scheme  notation.  Some  remarks.  The  remarks  are  con¬ 
cerned  mainly  with  our  treatment  of  the  improper  integral. 

Section  6:  Statement  of  difference  equations.  The  difference  equations  are 
set  up  under  assumptions  deemed  appropriate.  They  are  shown  to  be 
recursive. 

Section  7:  Outline  of  convergence  proof.  Estimates  deferred  to  Sections  8 
and  10  prove  that,  as  the  mesh  widths  vary,  the  solutions  of  the 
difference  equations  are  compact.  This  fact  and  uniqueness  (Theorem  2 
and  Section  2)  prove  the  solutions  of  the  difference  equations  to 
converge  to  a  weak  solution  of  the  corresponding  problem. 

Section  8:  A  bound  for  the  solution  of  the  difference  equations  and  an 

estimate  for  its  t-difference  quotients.  Two  of  the  estimates  applied 
in  Section  7  are  developed. 

Section  9:  Boundary  behavior.  The  manner  of  the  assumption  of  boundary  data 
and  related  questions  are  discussed  for  solutions  of  difference 
equations;  the  results  are  applied  in  Sections  10  and  11. 

Section  10:  Lipschitz  conditions  in  truncated  problems.  The  remaining 
estimates  required  in  Section  7  are  given. 
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Section  11:  Holder  conditions  in  untruncated  problems.  Proof  of  Theorem  3. 
Estimates  are  obtained  that  have  the  effect,  under  appropriate 
hypotheses,  of  showing  that  a  weak  solution  is  strong. 

2.  Uniqueness  and  continuous  dependence.  Uniqueness  and  some  allied 
properties  are  demonstrable  under  fewer  restrictions  as  to  K  than  those 
of  section  1.  Assumption  (iii),  namely,  can  be  replaced  by 

(iii)o  K(x,y,t,y')  is  integrable  with  respect  to  x,t,y'  on  S,p  . 
Furthermore, 

|K(x,y,t,y')[  <  Ko(y,y')  , 

where  Ko(y,y* )  is  integrable  with  respect  to  y'  ,  and  integrable  with 
respect  to  y  on  any  sphere  |y|  <  const.,  and  where 

^ KQ(y,y' )dy'  <  kQ  (kQ  =  constant)  . 

Our  earlier  uniqueness  statement  (Theorem  2)  is  included  in 

THEOREM  2.1.  Under  hypotheses  (i)  and  (iii)o,  a  function  u(x,y,t) 
belonging  to  V'(S,j,)  and  satisfying  inequality  1.7  with  f)  =  0,  g  =  0 
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is  zero  almost  everywhere  in  S^. 

This  fact  is  an  obvious  corollary  of  an  estimate  below  for  the  growth  of 


U(t)  =  sup  J(u(x,y,t)  )2dx  . 


In  stating  this  estimate,  for  convenience  we  suppose 

( i)Q  c  is  a  non-negative,  bounded  member  of  V" (S^) .  (Non-negativity 
is  merely  a  normalization  arising,  for  instance,  as  a  result  of  a  substitu¬ 
tion  u  =  e^*V  with  sufficiently  large  \.) 

THEOREM  2.2  (Continuous  dependence).  Under  hypotheses  (i)o,  (ii)> 
(iii)0»  (iv),  a  function  u(x,y,t)  belonging  to  V'(ST)  that  satisfies 
inequality  1.7  for  all  y  and  for  0  <  t  <  T  also  satisfies  the  integral 
inequality 

2.1  U(t)  <  (||(|>||'2  +  ||g||"2)e(2ko+  2)t 

for  0  <  t  <  T.  (The  norms  have  been  defined  in  Section  1.) 

Proof:  From  1.7,  since  c  >  0,  we  have 

2.2  (u(x,y,x) )2dx  <  ||4||'2  +  2  j'  ju^Kudy*  +  gu|  dxdt 

t  * 

for  all  vectors  y  and  all  t  in  the  interval  (0,T).  By  Schwarz*  in¬ 
equality  and  (iii)Q, 


(jKudy*^  <  (jlKldy*^  ^  J|K  |u2dy'^ 

<  k0 JK0(y»y')  (u(x,y',t))2dy'  ; 


hence,  by  Fubini's  theorem. 
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.3  J'  ^  J' Kudy'^  dx  <  kQ  j  dx  jKo(y,y* )  (u(x,y' ,t) )2  dy' 

^  =  kQ  fKo(y,y' )dy'  J (u(x,y!  ,t)  )2  dx 


<  k^U(t) . 


Consequently, 


2  ^  * 

J*  u  ^  $Kudy'^  dxdt|  <  ^  Ju2dxdt  J"  ^  J*  Kudy'^  dxdt  | 


<  k q  J* U(t)dt  . 

0 


This,  inequality  2.2  ,  and  Sehwarz'  inequality  applied  to  J* gudxdt  enable 
us  to  deduce  from  relation  2.2  the  inequality 

U(x)  <  ||<HI'2  +  2||g||"  ^  JVtJdt^  /2+  2ko  J  U(t)dt 

X 

<  I  I4>l  |f2  +  llgfl"2  +  (2k  +  l)  fu(t)dt  , 


valid  for  0  <  x  <  T  .  It  is  easily  seen  that,  for  0  <  t  <  T  ,  U(t)  < 
V(t),  where  V(t)  is  determined  by  the  conditions 


V'(t)  =  (2ko  +  l)V(t) 

v(0)  =  IH>||'2.  ||g||»2  . 


Thus,  estimate  2.1  is  justified. 
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3.  positivity,  Monotonic  dependence  of  solution  on  data,  coefficient, 

kernel,  and  inhomogeneous  part  of  equation.  When  the  coefficient,  the  in¬ 
homogeneous  part,  and  the  data  are  positive,  the  solution  will  be  positive* 
too.  This  is  more  precisely  stated  in 

THEOREM  3.1.  Under  hypotheses  (i)-(iv),  let  u(x,y,t)  denote  the  weak 
solution  belonging  to  class  W  of  which  we  are  assured  in  Theorems  1  and 


la.  If 


c  >  0,  g  >  0  ,  <J>>  0  , 


u  >  0 

at  almost  all  points  of  S^,. 

A  significant  consequence  of  positivity  is  that  u  depends  monotonical— 
ly  on  -c,  g,  K,  and  <j>  ,  as  we  assert  in 

THEOREM  3.2.  Consider  two  problems  of  the  form  specified  in  1.1,  1.2, 
and  1.3,  each  satisfying  hypotheses  (i)-(iv).  Distinguishing  corresponding 
quantities  in  the  two  problems  by  the  subscript  1  or  2  ,  assume  c^  >  0  and 

c2  -  °1  ’  gl  -  g2  *  K1  >  K2  ,  >  <j>2  ; 

denote  by  and  u2  the  weak  solution^  in  W,  of  these  respective  problems 
as  provided  by  Theorem  1.  At  almost  all  points  of  ST  , 


U1  i  u2 


Theorem  3.2  is  proved  by  noting  that  the  difference  v  =  u^-u2 
initially  is  non-negative  and  is  a  weak  solution  of  the  equation 
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Vt  + 


v 


=  jK^dy'  +  ((c.,  - 


°1>U2 


J<KrK2)u2dy'  +  «i  -  h) 


The  expression  in  parentheses  being  non-negative,  Theorem  3.1  applies. 

Theorem  3.1  need  be  proved  but  for  bounded,  Lipschitz-continuous  solutions 
since,  as  appears  from  Sections  2,  4,  and  7,  the  (weak)  solution  of  an  ar¬ 
bitrary  problem  can  be  approximated  by  solutions,  which  are  bounded  and 
Lipschitz-continuous,  of  related, approximating,  "truncated”  problems.  We 
shall  therefore  assume  that  u  is  bounded  and  Lipschitz-continuous . 
(Lipschita-continuity  with  respect  to  x,t  for  almost  all  y,  as  is  the  case 
for  a  strong  solution  in  the  absence  of  truncation,  would  work  equally  well.) 
Relation  1.4  holds  for  any  Lipschitz-continuous  solution,  and  we  may  make 
the  following  substitutions: 


F(u)  =  u  ,  u  = 


2k_6t 
e  0  v. 


-2kn6t 
<j>  =  e  ° 


where  6  is  an  arbitrary  number  >  1.  Thereby,  1.4  becomes 


3.1 


J  v(vt  +ZyrVr)dS 
D 


D 


ge 


— 2k  6t 
o 


} 


dxdt 


a  relation  for  the  new  dependent  variable  v.  Obviously,  it  suffices  to 


prove 


3.2  v  >  0  in  S. 


Imagine  that,  to  the  contrary,  for  some  positive  value- of  T  we  have 


3.3  m  =  inf  v(x,y,t)  <  0  . 

x  ejg 
all  y 
CKt<T 
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In  this  case,  a  vector  y*  exists  such  that 

3.4  m*  =  min  v(x,y*,t)  <  m/6  . 

x  e£ 

CKt<T 

Let  the  minimum  m*  be  realized  at  a  point  P*.  (x*,t*),  0  <  t*  <  T,  i.e., 
¥ 

3.5  v(x*,y*,t*)  =  m*  , 

and  assume  t*  the  least  value  of  t  at  which  this  minimum  occurs.  Thus, 


3.6  v(x,y*,t)  >  m*  for  0  <  t  <  t#  ,  x  t%  . 

By  relations  3.3  and  3.4  , 

jK(x,y,t,y')  v(x,y',t)dy'  >  k0m  >  k06  m*  • 

Hence,  and  because  g  >  0,  equation  3.1  implies 

cvdxdt  —  kQ6  I (2v  —  m*)dxdt  . 

D 


J*v(vt  +  Xyrvr)dS  >  -  J1  cvdxdt  -  kQ6  f( 


D 


If  N  is  a  neighborhood  of  P*  such  that  v(x,y*,t)  <  m*/2  in  n  N, 
we  thus  have,  since  c  >  0  , 


-3.7  jv(Vt  yrvr)dS  >  0 

S 

for  D  C  D  N  .  This,  however,  is  incompatible  with  the  minimality  of 
m*  expressed  in  3.5  and  3.6  ,  which  we  see  as  follows. 

,  consider  the  truncated  (d+l)-dimensional 

|yj|(t0  -  t.)  .  J  =  1 . a, 

t  <  t« 


For  t  >  t*  and  h  >  0 
0 


pyramid 


Pyr 


to,h‘ 


|Xj-x*|  < 

t*  -  3h  * 
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Knowing  P#  ,  we  can  select  h  such  that,  for  any  selection  of  P^ 
and  all  sufficiently  small  e  , 

3.11  Q  CR  . 

e,t] 

If  x*  is  interior  to  X  ,  this  is  evident  from  3.9  ,  for  h  in  this  case 

can  be  so  reduced  that  the  space  component  x  of  any  point  (x, t)  of 

Pyr^  k  also  will  be  interior  to  £ .  If,  on  the  other  hand,  x*  is  on  the 
o,n 

boundary  of  X  >  ve  see  from  1.3  and  3*5  that,  at  x*  ,  the  vector  y* 

points  outward  from  X  •  Hence,  the  characteristic - the  straight  line 

with  slopes  dx^/dt  =  y*,  i  =  1,  ...,  d  -  through  P*,  positively 

oriented  in  the  sense  of  increasing  t  ,  points  outward  from  R  s  the  seg¬ 
ment  of  the  characteristic  for  which  t*-h  <  t  <  t*  ,  in  particular,  is 
contained  within  R  .  Corresponding  segments  of  neighboring  (parallel) 

characteristics  also  will  be  contained  within  R  .  Q  ,  however,  is 

e»T) 

generated  by  such  segments.  Hence,  being  selected,  Q  will  be 

e,  t] 

contained  within  R,  as  asserted,  if  e  is  sufficiently  small. 
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Because  of  3.11  and  3.8  ,  we  can  choose  D  =  Q  in  3.7  .  In 

eJ  *1 

view  of  3.10,  inequality  3.7  then  reduces  to 

3.12  J  vdx  -  J  vdx  >  0  . 

F  F 

2  1 

Dividing  this  by  and  letting  e  0  proves 
v(P^)  -  v(p^**)  >  0  , 

where  P^**  is  the  point  with  the  coordinates  -  y*h  ,  t^  -  h  . 

1  i 

Now  letting  q  «♦  0  proves 


v(p*)  -  v(P**)  >  0  , 

where  P**  is  a  point  on  the  plane  t  =  t*  -  h  .  This  contradicts  the  mini¬ 
mality  of  m*  ,  i.e.,  contradicts  3.5  .and  3.6  ,  and  in  view  of  the  boundary 
condition  we  conclude  that  m*  =  0  .  This  means  statement  3.2  holds,  all 
we  needed  to  prove. 


4.  Truncated  problems.  Proof  of  Theorems  1  and  la.  A  problem  of  the 
type  described  in  section  1  will  be  called  "truncated"  if  g(x,y,t),  like 
K(x,y, t,y' )  ,  is  non-negative  and  both  are  zero  when  y  is  in  a  neighbor¬ 
hood  of  the  union  U  of  the  planes  y_  =  0  ,  s  =  1,  . ..,  d  .  Truncated 

S 

problems  play  an  essential  role  in  our  proof  of  Theorems  1  and  la,  for  their 
solutions,  properly  combined,  approximate  the  solution  of  any  problem,  while 
they  themselves  are  solvable  by  finite  difference  methods.  Their  solvability 
will  be  discussed  in  section  7,  their  approximative  capability  here:  the 
effect  of  the  combined  discussion  is  to  prove  Theorems  1  and  la. 
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Our  present  aim,  stated  more  precisely,  is  to  show  that  the  weak 
solution  of  an  arbitrary  problem  can  be  approximated  in  V' (S^,)  by  the 
solutions  of  appropriate  truncated  problems  with,  say,  continuously  differen¬ 
tiable  coefficient,  kernel,  data  and  inhomogeneous  part.  We  anticipate  the 
fact,  to  be  proved  in  Section  7,  that  such  truncated  problems  have  Lipsehita- 
Continuous  strong  solutions.  Three  steps  are  involved. 

The  first  step  is  to  understand  that,  under  any  stipulated  boundary  and 
initial  conditions,  equation  1.1  is  solvable  for  any  g  satisfying 
hypothesis  (ii)  if  solvable  for  any  such  non-negative  g  .  Any  g  subject 
to  (ii)  can  be  represented,  in  fact,  as  a  difference 

g  =  gx  ~  g2  (gx»  g2  - 

of  non-negative  functions  also  subject  to  (ii).  If  u^  and  u^  are  weak 
solutions  of  the  respective  relations 

=  gx  +^Ku1dy'  ,  Lu2  =  g2  +  J  Ku^y'  , 

where  Lv  =  v^  +  2y  v^  +  cv  ,  then  u  =  u^  —  u^  evidently  is  a  weak 
r  r  xr 

solution  of  the  equation 

Lu  =  g  +  J’Kudy1  . 

Thus,  it  suffices  to  be  able  to  solve  1.1  for  non-negative  g  ,  as  asserted. 
In  what  follows,  we  shall  generally  assume  g  >  0  . 

The  second  step  is  concerned  with  problem  1.1,  1.2,  1.3  when  g  >  0 
and  e,  g,  K,  and  ft  are  all  continuously  differentiable  besides  satisfying 
hypotheses  (i)  to  (iv);  for  convenience,  (i)Q  (Section  2)  also  is  assumed. 

The  (weak)  solution  of  such  a  problem  will  be  exhibited  as  the  limit  of 
solutions  of  appropriate  truncated  problems.  It  will  belong  to  W  .  More, 
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if  u  and  v  are  the  solutions  of  two  such  problems,  say  with 

(ut  +  Syrux  +  CjU  =  gx  +  jl^udy' 
r 

lu(x,y,o)  =  ^(x.y)  * 

|vt  +  *r\  *  V  =  g2+  iK2vd!r' 

(v(x,y,0)  =  £2(x,y)  , 

the  difference  u-v  will  be  shown  to  satisfy  relation  1.7  with  c  =  c^  , 

K  =  Kx »  g  =  gx  -g2  +  (°2  -cx)»  +  j^Ki“K2^v  dy'  »  f*  =  5  ** 

Theorem  2.2,  it  will  then  Bellow  that 

4.1  ||u-v||'2  <  C(|H)1-^2||*2+  ||grg2H"2+  ll(e2-c1HI"2  + 

+(max  C| —  K2 1 dy • ) 1 1 V 1 1 » ? ) 

sT  J 

where  C  is  a  constant.  (An  estimation  is  used  like  that  for  2.3.)  This 
step,  the  main  task  of  the  present  section,  will  be  carried  out  presently. 

The  third  step  is  to  perceive  that  the  solution  of  a  general  problem, 
one  satisfying  hypotheses  ( i j  to  (iv),  can  be  approximated  by  solutions  of 
problems  of  the  type  considered  in  the  second  step.  The  latter  solutions 
belonging  to  W  ,  their  limit,  the  solution  of  the  given  general  problem, 

too  belongs  to  W  :  Theorems  1  and  la  are  thus  proved.  We  shall  carry  out 

the  third  step  now,  assuming  the  second  to  be  valid,  and  then  return  to 
justify  the  latter.  We  make  assumption  (i)Q  for  convenience,  in  place  of 
(i),  and  take  g  >  0  ,  as  from  the  first  step  we  may.  Then  we  approximate 
{)  in  V’ (SQ)  by  a  sequence  of  continuously  differentiable  functions 
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k  =  1>2,  •  ••>  g  in  V"(ST)  by  a  sequence  of  non-negative, 
continuously  differentiable  functions  g^(x,y,t) ,  k=  1,2, ...  ,  and  c 
in  V"(ST)  by  a  bounded  sequence  of  non-negative,  continuously  differentiable 
functions  ck(x,y,t),  k  =  1,2,...  }  K  we  approximate  by  the  of 
hypothesis  (iiib).  It  follows  from  step  two  that,  for  each  k  ,  the  problem 
corresponding  to  {)^,  c^,  g^,  has  a  solution  u^  with  uniformly  bounded 
norm  j | u^ | | '  and  that,  for  all  k  and  m, 


+  (ek  +  Em)|lum|l'2)  ' 


Each  term  on  the  right  tends  to  zero,  however,  as  k,m-*  00 ,  because  of 
(iiib),  other  stipulated  approximation  properties,  and  the  existence  of 

the  u^  converge  in 
V' (ST),  and,  as  is  easily  seen,  the  limit  u  is  a  weak  solution  of  the 
given  problem  belonging  to  W  .  Theorems  1  and  la  thus  hang  on  the  validity 
of  the  second  step. 

In  justifying  the  second  step,  we  shall  actually  do  more.  By  V+ ( S^, ) 
will  be  meant  the  class  of  functions  w  ,  defined  on  S^,  ,  each  of  which 
is  the  limit  of  a  monotonically  increasing  sequence  of  functions  belonging 
to  V(ST)  .  Analogously,  by  V+(SQ)  will  be  meant  the  class  of  functions 
defined  in  Sq  each  of  which  is  the  limit  of  a  monotonically  increasing 
sequence  of  functions  belonging  to  V(SQ)  .  We  replace  the  earlier 
hypotheses  by  the  following: 


uniform  bounds  for  the  c,  and  the 

k 


II\H’ 


Hence, 
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(i)  c(x,y,t)  is  bounded,  and  -c  belongs  to  V  (S,p)  . 

(ii)+  g(x,y,t)  is  non-negative  and  belongs  to  V+(ST)  , 

(iii)+  (a)  K  >  0  . 

(b)  A  non-decreasing  sequence  of  non-negative,  continuously 
differentiable  functions  Km(x,y,t,y' ),  m  =  1,2,,..,  exists  such  that 

K(x,y,t,y')  =  lim  K  (x,y,t,y* ) 
m— *  oo 

at  all  points  of  . 

(c)  the  same  as  (iiic)  . 

(iv)+  $)  belongs  to  V+(Sq)  . 

These  assumptions  obviously  are  satisfied  when  g  >  0,  ({),  c,  g,  and  K  are 

continuously  differentiable,  and  (iiia,c)  hold.  Hence,  the  second  step  in 
the  argument  above  will  be  justified  when  we  prove 

THEOREM  4.1.  Under  hypotheses  (i)+  -  (iv)+  ,  a  weak  solution  of  class 
W  of  1.1,  1.2,  1.3  is  obtainable  as  the  monotone  limit  of  solutions  of 
appropriate  truncated  problems.  If  u  and  v  are  the  solutions  of  differ¬ 
ent  problems,  each  satisfying  hypotheses  (i)+  -  (iv)+  ,  u-v  is  subject  to 
an  inequality  of  the  form  4.1  . 

The  last  statement  is  an  obvious  corollary  of  the  first,  the  inequality 
in  question  being  obtainable  as  the  limit  of  analogous  inequalities  for  the 
(Lipschita-continuous)  solutions  of  the  approximating  truncated  problems. 

The  first  statement  depends  on  a  suitable  procedure  of  setting  up  truncated 
problems  to  approximate  the  given  one. 

Let  <p(s)  be  a  mono tonic,  infinitely  differentiable  function  on  the 
'positive  semi-axis  s  >  0  such  that 
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<p(s)  =  0  for  0  <  s  <  1/2 

=  1  for  s  >  1  . 

For  m  >  1  ,  the  function 

®  (s)  -  o(2s^m) 

m 

thus  is  monotonic  and  infinitely  differentiable  and  satisfies  the  conditions 

®  (s)  =  0  for  s  <  (l/4)m 

m  — 

=  1  for  s  >  (l/2)m  . 

Moreover,  for  s  >  0 

0_(  s)  /  1  as  m  /  co  . 

in 

If  p(y)  is  a  non-negative  function  of  the  vector  y  =  (y^,  ...,yd)  ,  we 
define  a  "truncation"  of  p  as 

d 

pm(y)  «  yp  ®m(  t )  •  p(y)  •  • 
i  =  i 

For  y  £  U,  a  truncation  evidently  tends  monotonely  to  its  "original": 

P„(y)  S  p(y)  as  a— tcD  if  y  |  D  . 

Let  us  now  consider  a  problem  satisfying  hypotheses  (i)+  —  (iv)+ 

We  shall  first  approximate  $)(x,y),  -c(x,y,t),  g(x,y,t),  and  K(x,y,t,y’) 
by  non-decreasing  sequences  {^(^y),  -c^(x,y,t),  gj£(x,y,t),  and 
K^(x,y,t,y' )  (k  =  1,2,...)  ,  respectively,  of  continuously  differentiable 
functions.  We  then  truncate  the  g^  and  with  respect  to  y  .  The 
truncated  functions,  denoted  by 
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gkm(x,y,t),  m=  1,2,  ..., 

Kkm^x,y,t,y, ^  m  = 

for  fixed  k  increase  monotonely  with  m  .  They  define,  for  each  k  and  m, 
what  we  have  called  a  truncated  problem.  Any  such  truncated  problem,  however, 
as  will  be  seen  below  (section  7),  can  be  solved  (by  finite  difference 
methods);  its  solution  u^m  is  bounded  uniformly  with  respect  to  k  and 
m  (Theorem  8.1),  is  Lipschita-continuous  and  thus  satisfies  appropriate 
relations  of  the  form  1.6  and  1.7,  and  increases  monotonely  with  m  and 
with  k  (Theorem  3.2) .  Hence,  the  limit 

u  =  lim  lim  u. 

,  km 

lr*  oo  m*oo 

exists  and  satisfies  1.6  and  1.7  ,  i.e.,  is  a  weak  solution  of  class  W 

of  the  original  problem.  This  proves  Theorem  4.1* 

The  remainder  of  the  paper  Will  be  concerned  with  suitably  smooth, 
truncated  problems. 

5.  Difference  scheme  notation.  Some  remarks.  Let  £  >  0,  fi  >  0, 
t  >  0  with  I  =  C A  an  integer,  and  set 
•6  =  t/4  . 

Let  L  be  a  positive  integral  multiple  of  t)  .  By  i  and  j  we  shall  mean 
multi— indices 

i  =  (i1,  ...»  ij),  -I  <  ig  <  I,  s  =  1,  .  v,  d, 

j  =  (jj»  l«j'sl  ^  s  =  1,  ...,  d, 

with  integral  components,  and  by  (  i  =  ((  i^,  ...,  (  i^)  and 
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r|  j  r]jd)  corresponding  points  of  and  of  y-space,  respectively. 

Symbols  such  as  i*  ,  £  i'  ,  j '  ,  r)  j  •  will  be  similarly  used.  For  any 
multi-index,  such  as  i  ,  we  shall  write 

j i I  =  ®ax  |ir|  . 

r 

The  set  of  all  multi-indices  i  such  that,  for  some  s  =  1*  ...,  d, 
i  =  I  ,  j i  |  <  I  for  r  +  s  , 

will  be  called  J)+  ,  and  the  set  such  that,  for  some  s  , 


1s  =  -1  ■  lipl  < 

I 

for  r  ^  s  , 

.  These  sets  correspond  to 

< 

and  B~  ,  respectively. 

ft  + 

Js 

<  0 

=  JF  for 

>  0 

=  empty  set  for  jg  =  0 


correspond  to  B  and  the  sets 

s,y 


a  _ 

d 

J 

-  U  ► 

^  r,j 

^  s 

r=l 

d 

j; 

-  v 

r?l 

^  r,j 

■  j:  *  j: 


for  jB  =  0 


to  B'  .  (Recall  that  u  is  prescribed  zero  on  B  and  is  unpree- 

9J 

cribed  on  B'  .) 

s,y 
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Let  S^rT',T  denote  the  set  of  lattice  points  (£i,  t)j,  nx)  in  xyt- 

space,  where  i  and  j  are  multi-indices  as  above,  and  n  is  a  non-negative 

integer.  Our  aim  in  the  difference  scheme  below  will  be  to  approximate  the 

solution  u(x,y,t)  of  a  given,  truncated  problem  by  functions  u^,r',T(x,y,t) 

L 

defined  on  ^,r*,T  .  When  5  ,  f)  ,  x  ,  and  L  are  regarded  as  fixed,  for 
brevity  we  write 

u£,T1,T(S  i,  r| j ,  nx)  =  u"j  , 

and  for  any  function  ^(x,y,t),  unless  otherwise  specified, 


'KU,  r)j,  nx)  =  t|>£.  ,  \|>x  Ui,  no ,  nx)  =  ^  ,ij  * 

s  s 

etc.  In  the  case  of  the  kernel  K  ,  we  write 


vn  _  — lr— d  -2d 
5  11 


K(x,y,t,y')  dy'dxdtdy  , 


Ti 


ij  v 


the  domain  of  integration  being  the  parallelepiped 


<  X  <  *(i  +  1)  ,  s  =  1,  ...»  d  , 

5  ™  S  S 

ni8  <  ys  <  n(js  +  i)  ,  " 

n 

^iii'  :  ni;  <  y;  <  n(i;  +  D  , 

nx  <  t  <  x(n  +  1)  . 


(In  this  definition,  it  is  presumed  that  K  is  defined  for  -I  <  x  <  I+E  . 

—  s  — 

We  might,  for  instance,  take  K  =  0  for  any  x  >  I . )  When  all  the  defini— 

s 

tions  of  this  section  are  concluded,  we  shall  come  back  to  note  certain 
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consequences  of  this  one. 

Let  Xr  be  the  translation  operator  whose  effect  on  the  multi-index  i 
is  to  increase  its  r-th  component  by  one  while  not  changing  the  other  com¬ 
ponents: 


Xr^  ~  W  •••*  S-l'  ir»  Vt-l*  V 


=  (i1,  ...»  ir  +  1,  ir+1»  •••>  id)  • 


The  effect  on  a  function  indexed  with  i  ,  say  ,  is  defined  by 


v  ,n  ,n 

Xr  +ij  ~  +Xji),J  * 


The  inverse  of  this  operator,  denoted  by  X~^  ,  reduces  if  by  one  while 
leaving  all  other  indices  the  same.  The  operator 


Tr,j=Xr  for  Jr<0 

=  X"1  for  j  >0 
r  r 

has  the  effect  of  moving  i  nearer  $  ,  . 

r»  J 

As  analogous  to  Xr>  let  Yr  denote  the  translation  operator  whose 
effect  is  to  increase  by  one  while  not  affecting  other  indices,  and 

Yr^  its  inverse.  In  terms  of  these,  conventional  *-  and  y— difference 
operators  sure  introduced  as 


5.1 


and 


5.1 


r  =  1, 


d  , 


. . ., 


d  , 


29 


NOLTR  62-193 


but  we  also  shall  use 


5.2 


and 


6  . 
r,J 


1  -  X"1  when  j  >  0  , 

r  r 

X  -  1  when  j  <  0  , 

r  "r 


5.2 


X  -  1  when  1  >  0  , 

r  r 

1  -  X"1  when  j  <  0  . 

r  r 


We  note  that 


5-3  Jr  5r>j  =  |jrl  (1  -  T  >  ; 


this  property  will  turn  out  to  be  a  significant  one. 

We  now  return, as  promised,  to  our  Jcernel  approximation.  Anticipating 
hypothesis  (iii)1  in  section  6,  assume  K(x,y,t,y' )  to  be  continuous  at 
almost  every  point  of  .  Note  at  the  same  time  that  the  step  function 

KL,fl’T(x,y,t,y')  =  in  for  |j|  ,  [j'f  <  L 

=  0  for  | j |  >  L  or  [j'|  >  L 

tends  to  K(x,y,t,y')  in  £T  at  each  of  its  points  of  continuity,  as 
5,  n»  T-*0,  L-*oo  .  Then 


5.4  lim  K^,T|,T  (x,y,t,y* )  =  K(x,y,t,y' )  almost  everywhere  in  2T_  . 

L  T 

L  oo 


Furthermore, 

-5.5  J4»,'»'t(x,y,t,y»)  dy«  <  jK(x,y,t,y» )  dy»  <  ko  . 
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These  facts  suggest  that  we  attempt  to  approximate  the  improper  integral  in 
equation  1.1  by  sums  (these  are  essentially  Riemann  sums)  of  the  form 


d 

*1 


I 

|j'i<L 


The  attempt  succeeds  when  the  step  functions 


4,ll,T(x,y,t)  = 


n 


for 


nr  £  t  <  (n  +  1)t 
(1s+  1)5 

1is<Is<  Dl 


(s=l,...,d) 


are  uniformly  bounded  and  tend  to  a  limit  u(x,y,t)  at  almost  all  points 
(x,y,t)  as  L->  cd  and  Z,  q,  t  -♦O,  the  approach  of  each  parameter  to  its 
limit  being  through  suitable  values.  To  be  more  specific,  let  q^, 

L^j  k  =  1,2,...  ,  be  sequences  of  values  of  5,  q,  x,  L,  respectively, 
such  that  the  functions 

^(x^t)  =  u^  X(x,y,t)  , 

as  stipulated,  are  uniformly  bounded  and  tend  to  u(x,y,t)  almost  everywhere 
as  k  oo  .  In  terms  of 

Kk(x»y>t,y' )  =  (x,y,t,y' )  , 

the  Riemann  sums  considered  can  be  written  as 
jKk(x,y,t,y')  u^x^tjy')  dy'  . 
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We  shall  need  to  prove  the  existence  of  the  limit  of  these  sums  as  k— *  oo 
only  in  a  weak  sense  to  be  described  below.  Let 

Qab:  as  <  ^8  ^  bs  ’  s  =  1»  d  ' 


be  a  parallelepiped  in  y-space.  Let  v(x,t)  be  a  continuous  function  which 
vanishes  for  x  and  for  t  >  T,  where  T  is  an  arbitrary  positive 

constant.  We  shall  need  to  know  only  that,  under  the  hypotheses  explained, 
for  any  V  as  described  and  any  numbers  a  and  b  ,  we  have 


5.6 


J'  dy  vdxdt  J”'  (K^u^  -  Ku)dy* 


*ab 


as  k  -#  00 . 

To  see  this,  we  write  the  integrand  as 
K(\  -  u)  +  Uj^  -  K) 

and  express  the  above  integral  as  the  sum  of  two,  the  first  correspond¬ 

ing  to  the  first  summand  above  and  the  second  2^  corresponding  to  the 
second  summand.  The  boundedness  of  jujJ ,  the  convergence  of 
and  5-5  show  that  J^k  and  both-*  0  as  k-*a>  .  Thus,  contention  5.6 
is  proved. 

We  also  remark  that,  from  the  definition  of  and  hypothesis  (iii), 

5-7  nd  I  <  ko 

y 

(this  is  the  same  statement  as  was  made  in  5.5  .)  When  assumption  (iii)^ 

below  is  made  concerning  difference  quotients  for  K,  the  following 
.-analogous  statement  can  be  made  about  , :  if  6  stands  for  6 , 


\  and  1^, 
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er/T)  ,  or  the  analogous  differencing  operator  with  respect  to  t  , 
5.8  1.^.1  <  ^  - 

V 


6.  Statement  of  difference  equations.  The  difference  equations  set  up 
below  will  be  justified  only  for  truncated  problems  in  which,  besides  assump¬ 
tions  (i)Q  ,  (ii),  (iii),  (iv),  (v)^  of  sections  1  and  2,  the  following 
additional  hypotheses  also  are  satisfied: 

(i)^  c(x,y,t)  and  g(x,y,t)  are  Lipschita-continuous  with  respect  to 
x,y,t  with  uniform  constants  we  shall  denote  by  nnH  ,  respectively. 
Furthermore,  0  <  c  <  cq  and  | g |  <  bQ  ,  where  cq  and  bQ  are  constants. 

K(x,y,t,y')  is  continuous  at  almost  every  point  of^^  •  If 
DK(x,y,t,y' )  symbolizes  a  difference  quotient  of  K(x,y,t,y')  with  respect 
to  any  of  the  2d  +  1  arguments  xr,  yg,  t,  r,s  a  1,  ...,  d,  e.g., 


*')  -  K(x, 
At 


(At  >  0)  , 


we  have  for  a  suitable  constant 


k 


1 


A 

|DK(x,y,t,y» )|dy'  <  . 

yJ 

(iv)^  $)(x,y)  is  Lipschita-continuous  with  respect  to  x  with  a 
uniform  constant  denoted  by  Also  |{>|  <  where  is  a  constant. 

(iv)  i>(x,y)  is  uniformly  Lipschita-continuous  with  respect  to  y  . 

We  now  describe  our  difference  scheme.  Corresponding  to  1.2  is  the 
initial  condition 
6-1  =  fa  i,qj)  , 
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and  corresponding  to  1.3  the  boundary  condition 


u,  .  =  0  on  •  i  s  —  I*  •••jd 

ij  s>  J 


Corresponding  to  the  integro-differential  equation  1.1  we  set  up  the 
difference  equations 

1/  n+1  n  w  V  •  1  r  „n+l  „n 
T(Uij  *  “ij1  +  V  1  h  X  6S,jUij  +  - 


n  d 
=  +  i 


\  vn  n 
4—  Kijj'Uij' 


j'l<L 


for  the  following  values  of  the  indices: 


n  =  0,1,  ...;  | j  j  <  L  ;  ji[  <  I  or  i  £  $  . 


This  scheme  is  recursive:  for  each  choice  of  n  and  j  ,  the  values 
n+1 

of  u„  ,  appropriately  ordered,  Eire  successively  obtainable  from  previously 

calculated  quantities.  This  fact  will  be  more  apparent  from  a  rewriting  of 

equations  6.3  in  terms  Of  the  T  .  .  We  simply  apply  5-3  in  6.3  . 

s,  j 

Then,  rearranging,  we  have,  as  equivalent  to  6.3, 


d  +  «  H  Ijgln)^1 

S  J 


(1  -  TcJjJugW 

s 


n  .  d  \  vn  n 
+  T%j  2-  Kij3.  “ij.  • 


llM<i 


Assuming  the  u’jjV  to  be  known  for  fixed  n  and  all  i  and  j  ,  we  shall 
see  from  this  how  the  u^  axe  determined.  The  determination  is  immediate 
in  the  case  j  =  0.  In  the  contrary  case,  j  f  0,  u*  is  expressed  in  6.5 
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in  terms  of  quantities  indexed  with  n  ,  presumed  known,  and  quantities  of 


the 


n+1 


type  u. , ,  for  which  i*  falls  one  step  nearer  some  .  than  i 
i  j  s,j 


(the  ^  . ,  we  recall,  carry  boundary  data).  To  be  more  exact,  consider  a 

s#  j 

particular  j  ^  0.  The  intersection 

A  'Vi 

Js*°  ,J 

of  all  the  ^  .  for  non-vanishing  j  is  not  empty;  let  i#  denote  any 

s>  3  s  J 

index  belonging  to  this  intersection,  and  set 

i'  =  (TV  T”1 . ) i*  , 

J  j  1  o  s>*  i 

Js  T 

the  product  in  the  right  member  to  include  every  operator  T  \  for  which 

s»  J 

Js  4  0.  As  i*  varies  over  all  possible  positions  and  kg,  s  =  1, ...,d, 
over  all  needful,  non-negative  integers,  the  index 

-k 


(  TT  T  B)V 
js*  o  i 


varies  over  all  the  lattice  points  i  at  which  is  not  already 

prescribed  in  advance.  Thus  it  suffices  to  show  that,  for  any  i*  ,  the 

J 

quantities 

6.6  (  TT 

j*  +  o  ,J  V 


can  be  recursively  determined.  First,  u?^  =  (  'll'  T~x.  )u?tt 

V  jo+0  V 


%  n+1 


can  be 


calculated  from  6.3  ,  being  known  on  .  when  j  ^  0  .  From 

s,j  s 

this  result,  >  for  ^  °»  similarly  can  be  found 

J 
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Induction  would  show  that  all  the  quantities  6.6  ,  in  suitable  order,  can 
be  determined  by  appropriately  continuing  this  procedure.  Hence,  the  scheme 
6.1-3  is  recursive,  as  asserted. 

For  future  reference,  we  note  that  equations  6.3  also  may  be  written 
as 

6.7 


'ij'Uij 


3  .A  U,  , 

sj  s  ij 


n 


n 


j'l<L 


KiJJ •  “ij  > 


where,  for  j  >  0  , 
s  *■ 


6.S.  Bsj  =  0,  Csj  =  , 


and,  for  j  <  0  , 
us 


6.8b  Bgj  =  .  csj  =  0. 


7.  Outline  of  convergence  proof.  The  proof  described  in  this  section 
applies  to  any  truncated  problem  satisfying  assumptions  (i)-(iv)  and  (v) 


(£) 


of  section  1  and  (i)^,  (iii)^,  (iv)^,  (iv)^  of  section  6.  For  such  a 
problem,  the  solutions  of  the  difference  equations,  u^,fl,T(x,y,t)|  are 

equibuunded  and,  for  any  £  >  0,  satisfy  Lipschitz  conditions  in  ST  -  S 

with  respect  to  x  ,  y  ,  t,  r,  s  =  l,...,d,  which  are  uniform  in  this  set 
r  s 

for  all  possible  x,y,t,£,q,x,L.  This  will  be  proved  in  sections  8  and  10 
below.  Because  of  the  uniform  bounds  and  the  uniform  Lipschitz  conditions, 
from  any  infinite  set  of  lengths  L  tending  to  infinity  and  of  suitable 
mesh  widths  K,  t|,  t  tending  to  zero,  sequences  k=l,2, ...» 
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can  be  selected  such  that  the  solutions 


u^(x,y,t)  =  u.k  k  k(x,y,t) 
Lk 


converge  continuously  ^  to  a  function  u(x,y,t)  that  is  Lipschitz  continuous, 

(  c  ) 

uniformly  so  in  S^,  —  S'  '  for  any  e  >  0.  This  limit  function  will  be 

seen  to  be  a  strong  solution  of  the  problem. 

By  continuity,  u  plainly  satisfies  the  requisite  boundary  and  initial 
conditions,  and  we  shall  now  see  that  equations  1.1  hold  almost  everywhere, 
as  well.  Let  v(x,t)  be  a  function  of  class  C'  in  the  region 

x  e  £  ,  0  <  t  <  T  (T  >  0) 

that  vanishes  identically  in  a  neighborhood  of  the  boundary  of  the  region, 
and  let  v?  =  v(£  i,nx)  denote  its  restriction  to  the  lattice  = 


VTk 


On  this  lattice,  consider  the  difference  equations  6.3  deter¬ 


mining  u^  .  We  multiply  the  members  of  equation  6.3  by  ,  sum 

over  n  and  i  ,  and  eliminate  u-differences  by  Abel's  method  to  obtain 


_  K  un  i  (vn+1  -  vn)  -  T~  n  j  u1*1  i  5  .vn+1  +  Y  cn.vf  V,  = 

Uij  Tk^Vi  V  L-  Vs  ij  s,j  1  i—  1J  1  1J 


n,i,s 


Zn+l  n  .  ,  n+1  d  Kn  n 

vi  ««  +  2_  vi  \  2_  Kijj'  v  • 


We  now  multiply  this  by  t  x.  t)  ,  for  arbitrary  a  and  b  sum  over  4 

k  K  k 

such  that  j  t)  g  Qa  b,  and  let  k->  oo  .  Continuous  convergence,  and  the 
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i 


remarks  in  section  2  concerning  the  improper  integral,  show  that,  as  the 
result, 


^  y  v  +  cv)dxdt 
s  s 


dy/ \  V^8  + ^ 


Kudy' )dxdt 


*a,b 


Thus,  u  satisfies  equation  1.1  in  a  weak,  integral  sense.  For  any  e>0, 

(  £  ) 

however,  u  is  Lipschitz  continuous  in  S  -  S  .  If  Q&  therefore, 

intersects  no  coordinate  plane  y  =  0,  the  partial  differentiations  in  the 

s 

above  formula  may  be  transferred  to  u  by  integrating  by  parts,  and  we  thus 
obtain 


a  formula  valid  for  every  fraction  v,  and  every  a  and  b  ,  as  described. 
It  follows  that  equation  1.1  holds  almost  everywhere,  and,  hence,  is  a 
i  strong  solution,  as  contended. 


8.  A  bound  for  the  solution  of  the  difference  equations  and  an  estimate 
for  its  t-difference  quotients.  In  this  section,  under  the  hypotheses,  say, 

of  Section  6,  we  shall  first  prove 

8.1  |u£,T1,T(x,y,t)|  <  M(t)  , 

kt  b  k  t 

where  M(t)  =^e°+^2(e°—  1)  .  Then  we  shall  obtain  a  Lipschitz 

o 

condition  with  respect  to  t  of  the  form 
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ejc^esgjwpwiBtaRpF'S' 


F?W*re«  Jffj  *■  ^«^5JHiqM*Sf3K,7! 


8.2  |t£»n,T(x,y,t‘)  -  4*T1»T(x,y>t«)|  <  A(1  +Z|ysl)  (t*  -  t"|  , 

s 

where  A  depends  on  T,M(T)»  and  the  constants  in  hypotheses  (i)^,  (iii)^, 
(iv)^,  (v)^  . 

For  0  <  nx  <  T  ,  set 
Mn  =  max|u^j  f  , 

the  maximum  being  taken  for  |  i  |  <  I  ,  |  j  |  <  L  .  From  the  difference 

equations  6.5*  the  estimate  5*7  ,  and  assumptions  (i)Q  and  (ii)  ,  we 
have 

<1-M>XIJ>>|i£1|  <  «„+  «ZUsl'>»n).l  +  V  +  tkoMn  • 

s  s 


I  1 

For  suitable  determinations  of  i  and  j  ,  however,  |u?T  f  =  M  .  . 

|  •  n+J- 
\  Hence,  we  arrive  at  the  inequality 

? 

r 

i 

i 

I  <  (1  +  k  t)M  +  b  x  , 

|  n+1  —  o  n  o  * 

from  which,  by  induction, 

i 

M  <  (1+k  T)nM  +  b  x[(l+k  x)""1  +  (1+k  x)“"2  +  ...  +  l] 
j  n  — •  o  o  o  *•  o  o  J 

t 

j 

}  k  nx  b  k  nx 

1  <  Mce°  +r(e  0  -i)  • 

o 

I 

Since  Mq  <  ^)q  ,  inequality  8.1  follows. 

To  prove  8.2,  we  shall  consider,  in  an  analogous  way,  difference 
-equations,  initial  values,  and  boundary  values  pertaining  to  the  t-difference 
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quotients 


n 

"ij 


.  n-t-1  n 

!u  ~uij 


(Difference  quotients  with  respect  to  x  or  y  cannot  be  estimated  in  quite 
the  same  way,  the  first  because  of  the  boundary  condition,  the  second 
because  the  pertinent  difference  equations  will  contain  x— difference  quotients 
not  known  to  be  bounded.)  Boundary  values  are  zero,  as  before.  Initial 
values  satisfy  the  inequality 


8.3 


I’ll  1  5  CJ  * 


where  C.  =  (j)E|j|»l  +  b  +  (c  +k)(b  ,  as  we  easily  see  after  sub- 

1  *  s  o  o  o  o 

s 

stituting  ujj  =  w^  +  u°^  in  equation  6.7  (n=0).  The  resulting  relations 

may  be  written,  in  fact,  as 


d+  «  E  I  j  U?.  +  «  E  B.U.  +  «  E  CH.X”1w. 


’s'^-ij 


'u  y 


•  A  "•  1 

SJ  s  ij 


sj  S  "ij 


+  ^  ^  Bsj(~Uij+XsUij)  +e  21  Csj  (~uij+xsluij  * 
s  s 


,o  d  Vr  vo  o 

+  T*U  ^  Ki3j'  UIJ  • 


Hence,  for  W.  =  max  |w. .(  ,  we  have 
J  ^  ij 


d  E  lJsl^i)|wijl  <  «  £  |jsh  Wj  +  Cjt  , 
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and  therefore 

w3  *  V  ' 

this  was  the  contention  in  8.3  . 

Having  considered  the  boundary  and  initial  values  of  v"^  ,  from 

equations  6.7  we  must  now  derive  the  difference  equations  the  latter 

quantities  satisfy.  These  difference  equations  are  easily  seen  to  be  again 

of  the  form  6.7  ,  their  coefficients,  and  the  kernel  in  the  Riemann  sum, 

being  in  fact  unchanged,  and  the  new  inhomogeneous  term  that  takes  the  place 

of  g1?,  consisting  of  bounded  quantities.  Dividing  by  1  +  C.  ,  we 
J  J 

obtain  a  system  of  relations  that  we  will  regard  as  equations  for 

v”  =  v!1,  /(I  +  C.);  by  (v).  (Section  1),  the  kernel 
ij  13r  J  4 


IiV1  +  cj'1/<1  +  cj) 

in  these  equations,  acted  on  by  r|d  T  ,  gives  a  uniformly  bounded  result. 

J1 

n 

The  boundary  values  of  the  are  zero,  their  initial  values,  by  8.3  , 

bounded.  Hence,  by  the  argument  by  which  we  proved  8.1  ,  is  uniformly 

bounded,  i.e.,  (1  +  -  u^-I/t  is  uniformly  bounded,  at  least 

for  0  <  m  <  T  .  This  implies  statement  8.2  . 
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The  results  do  not  depend  on  any  kind  of  truncation.  We  also  shall  estimate 
the  boundary  values  of  x-difference  quotients  of  the  solutions,  but  these 
only  in  truncated  problems.  Three  theorems  are  formulated,  based  on  the 
hypotheses  of  Section  6. 

THEOREM  9.1.  If  yr  +  0,  we  have 

9.1  |i4>T1’T(x,y,t)|  <  Ar>y|xr  -  Cr>y|  for  0<t<T  , 


where  A  =  maxC^  ,  (b  +  k  M(T))/|y  |).  (C  is  defined  in  Section  1, 
r,y  loo  r  r*y 

M(T)  in  Section  8.) 

Proof:  Select  any  fixed  r  among  the  indices  from  1  to  d  .  With 
fixed  T  >  0,  for  jr  |  0  set 

P.  =  P  .  =  max  |T~m.u?.| 

jm  rjm  i£^  1  r,j  ij 

J 

<Xn<T/r 

From  equations  6.5  ,  for  (n-t-l)x  <  T  we  obviously  have 

9.2  <!♦«  I  |JSI-|)|T^1|  <  +  «  r  U.InPj. 

s^r 


+  ■©  j  fi  P.  ,  +  br  , 
Nr'  1  j,m-l  * 


where  b  =  b  +  k  M(T)  . 
o  o 

Either  |  reaches  its  maximum  P.  initially*  in  which  case 

r»j  ij  Jm 


9'3  V=^  .|Tr,jUiJl  *  *1*  ’ 

r»  j 
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» 


provided  9.4  holds  for  m"  >  m  >  m'  .  In  the  last  inequality,  replace 
m"  by  m,  and  let  m'  be  the  least  possible  integer.  If  m'  =  0,  we  have 


V  5  W\ir\<\  » 

since  Q  =  0  from  the  boundary  condition;  this  implies  9.1.  If  m'  >  0, 
relation  9.3  must  hold  for  m  =  m*  :  hence. 


PJ-  *  V  5  +lljn 

<  m£*  max  (^  ,  t/|jrh); 


this  again  implies  9.1. 

THEOREM  9.2.  Let  0  <  a  <  1,  and  let  e  be  a  positive  number  divisible 
by  r)  such  that 
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9.5  e  ^  <  b0  +  kQM(T)  . 

Then 

£/n  r 

9.6+  f)  J  max  |ujj,^,T(x>y^»  •  •  •»yp_2*^rT',^r+l*  *  '  — 

4-=-,  (x,t)etT 
Jr 

<  C  (C  +  x  )Q 
—  a  r 

and 

9.6  ri  y  max  ’ ’^r-L^r^’^r+l*  *  ’  *  — 

Hi  <x-t>,JET 

Jr 

4Ca(C  -  xr)a  , 

where 

C  sO-^naxMD.J&C)  +  (b+k  M(T))  [(2C)1“a(l+log(l+ee)) 
a  j.  o  u 

+  max  C1  alog(l  +  C  a)] 

CKC<2C 

Proof:  The  following  calculations  presume  0  <  |jrh  <  e  .  Let 

Qjm  =  Bax  {Pjm»  V  °  ‘ 

We  note  that 

«|jrh(4)1m^)  <  «|jrh(<>1(®-l)^)  +  b  t  , 
since  -  «|jrh  ^  5  +  br  =  (-|jrh  ^  +  b)-t  >  (b -e  >  0 

by  assumption  9.5;  hence, 

9.7'  «|jrh(^)  <  «lJrh  Qj>in_i  +  b  x  • 
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In  the  case  in  which  9.3  holds,  therefore,  we  obviously  have 
9.7"  «|JrlnPJm<  «IJrll«3,m.1  +  b7  • 

Even  when  9.3  fails,  since  9.4  then  is  valid,  inequality  9.7"  again 
can  be  deduced.  Thus,  9.7"  holds  without  exception.  With  9.7'  ,  it 
proves 

™  <  «UrllQj^.1  7  b,. 

We  shall  rewrite  this  result  as 

(1  +«|JrU>Qjm  <  • 

Setting  aj  =  dr  j  =  (1  +  ■©Ijj.lq)”1  ,  we  thus  have 

9'8  V  5  aJQJ»  +  (1  -  +  “jb7  • 

Now  define 

e/r] 

V  =  "  21  • 

if1 

To  prove  9.6  ,  it  suffices  to  estimate  AQm  appropriately.  Multiplying 
9.8  by  T|  aj  ,  p  =  0,1,...,  and  summing  over  jr  gives  us 

Aom  ^  Alo  +  A0,ro-l  -  Ai,m-1  +  ba°*(1  +  *  e) 

Apm  <  Vl,m  +  Ap,m-1  "  Vl,m-1  +  »  p  =  1»2»***» 

since 


45 


NOLTR  62-193 


e/t) 


e-j 


p+i 


dy 


Jr=1 


(l+Sy) 


|  log  (1  +  -Qe)  for  p  =  0 


<  pS  for  P=  1»2 . 


We  rewrite  the  last  results  as 


9-9o  Aom  -  w  S  Al»  -  Al,m-1  +  *  1oa  (1  +  «*> 


and 


9,9p  APm  ~  Ap,m-1  ^  Ap+l»m  “  Ap*l,»-1  +  **/p  ’  p  =  1»2’ 


Noting  that  Aqq  =  A  q  =  0  ,  p  =  1,2,...,  by  summing  with  respect  to  m 
we  have 


/'I  .... 


9.10  A  <  A.  +  bn*  log  (l+0e) 
o  om  —  lm  p 


and 


9,10p  Apm  <  Ap*-l,m  +  *  p~  1»2»“*  * 


Hence, 

9-u  Aom  <  Vl,m  + « (l0* (Me)  ♦ 1  ♦  I  ♦  -  ♦  0 


We  now  note  the  estimate 


Ap+l,m  <  « 


f  — <  Qto 

J  (l^y)^1 


(this  assumes  p  >  l),  where  Q  =  max  (M(T),2^C)  .  Also,  we  recall  the 
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fact,  known  from  Euler,  that  the  expressions 
1  *  3  *  •  *  *  +  p  ~  ^ 

are  positive  and  less  than  1  .  From  9.11  we  thus  have 

9.12  A  <  Q/9p  +  boS  (log  p  +  log  (l  +  -fie)  +  1)  . 

om 

Now  choose  p  such  that 

“  <  (■*)*  <  if  0  <  m£  <  1  , 

p  =  1  if  m?  >  1  . 

We  thus  also  have  that 

log  p  <  log  (1  +  (m K)~a)  i 
hence,  inequality  9.12  leads  to 

*om  <  !<"*>“  +  ba<«>°  - 

where  b  =  b(2C)1_a(l+log  (1+ee))  +  b  max  H1*®  log  (l  *  tf°)  . 

°  CKC<2C 

This  result  implies  9.6. 

THEOREM  9.3.  For  a  particular  index  r  ,  suppose  K(x,y,t,y')  and 
g(x,y,t)  to  be  zero  when  |yr!  <  u  ,  u  being,  say,  an  exact  positive  multiple 
of  tj  .  Then  for  Jr  4  0  ,  i  e  ^  ^  ,  0  <  qx  <  T  , 
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(t;  1  -  i)uu 


where  “  “ax  ,  <bQ  +  koM(T))/e,)  . 


Proof:  Take  0<  |jr|fi<u  and  m=l  in  9.2.  Since  P 


Pj,0  =  0  »  and»  under  present  hypotheses,  b  can  be  replaced  by  zero, 


obtain 


9iU  (1+4>£  Usll)H;^1|  <  U*«  £  |j  1,)F  . 


Either  jT^.u^l  reaches  its  maximum  P  initially,  in  which 


?Jl“H  |Tr>ijl  -<  *1.<  > 


or,  for  so.,  value  of  n  ,  |t;\|uJ*1|  an  assuae  thle  maxim™.  .  In  the 

letter  eaee,  w.  can  replace  la  9.14  hr  P,,  ,  aad  the  reeolt 

|Tr,jUij  1  J 

will  prove  PjX  =  0  .  Therefore,  P^  <  ^  without  qualification  when 
0  <  IJrh  <  w  ,  while  Theorem  9.1  applies  in  the  contrary  case  |j  (t,  >  u  . 
Thus,  9.13  is  completely  proved. 

10 *  Mpschitz  conditions  in  truncated  problems.  Our  proof  of  convergence 
for  truncated  problems  (section  7),  under  assumptions  (i-iv),  (v)  ,  (i) 

(iii)^,  (iv)1,  (iv)y  (sections  1  and  6),  presumed  knowledge  of  Lipschitz 

.conditions  with  respect  to  xr,  y  ,  t,  r,s  =  1,  d,  satisfied  for 
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0  <  t  <  T  by  u^’^’^Cxjyjt)  uniformly  in  x,y,t,5,ri,x,L.  A  suitable 
Lipschitz  condition  with  respect  to  t  already  is  to  be  had  from  section  8. 
Here,  we  shall  deduce  the  other  Lipschitz  conditions  needed,  proving,  namely, 
that,  for  0  <  nx  <  T  ,  the  difference  quotients 

CVJjl  - 

are  bounded  uniformly  with  respect  to  (,  q,  x,  L  and  all  indices. 

Let  u  describe  the  truncation:  i.e.,  suppose  K(x,y,t,y' )  and 
&(x,y,t)  to  be  zero  if,  for  any  s,  |y  |  <  u  ;  u  >  0  .  First  considering 
differences  with  respect  to  x  ,  i.e.,  i-differences,  we  shall  prove 


i  „  k  nx 

S-1|5ruijl  <  e  0  max 


b1+M(nx)(c-+k, )  k  nx 
-f  1  1  1  (e  0  -  1)  , 


being  the  constant  defined  in  Theorem  9.3. 

Estimate  10.1  will  be  obtained  from  difference  equations  for  the  i— 
differences  considered,  which  we  may  obtain  by  applying  6r,  defined  in  5.1, 
to  both  members  of  equations  6.7  .  We  use  the  identity 


6  (a.v. )  =  a.  6  v.  +  X  v. 
r  i  i  i  r  i  r  i 


valid  for  any  quantities  a^  and  v^  indexed  with  i,  to  obtain 
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u  +  ij.uir1  vI}1  =  a  -  «y>r*  6ru; 


H  \  t 


*-•  %  Vs  Chr •<?  +  *  S  V?  ChA?  +  T  5'V?J 


.  d  r-  vn  r-l,  ,  n  .  -n 
+  T  1  KijJ-  5  sruy  +T  S1J  > 


where 


S?J  —  rl(6rC«)ZAj  +  rl<»rKUi')XAj' 

By  hypothesis  (i)^  and  inequalities  5.8  and  8.1  , 

10.3  |S“  |  <  (Cl  +  kjMw)  . 


Because  of  the  existence  of  the  boundary,  there  are  i—  differences 

6  u^T1  =  ufT1  ~  X  u^1  in  right  members  of  equations  10.2  that  never 
r  ij  ij  r  ij 

appear  on  the  left.  These  are  of  two  types,  the  first  being  of  those  with 
indices  i  e  J  ,  ,  s  4  r;  these  vanish.  The  second  type  is  of  differences 

s»  j 

for  which  i  or  X  i  is  on  ^  t  •  We  shall  speak  of  the  differences 
r  r,  j 

6  u^t1  of  this  type,  or  of  the  corresponding  difference  quotients,  as  being 
r  ij 

"situated  on",  or  "on",  ^  .  .  Their  behavior  is  controlled  by  inequality 

r»  J 

9.13,  which  we  may  reword  as 
10.4  f’Vyl  < 


For  0  <  nx  <  T,  set 
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Bi  =  -?  «_1l6r“yl  •  Dn  =  *“  (Di>  V  • 


From  10.2  and  10.3  we  have 


io.5  a  +  «  2  |jaln)s"’1l®ruJJ1l  <  u  +  v)d;  +  «srijsl,iDA+i+b,T 


where  b'  =  b^  +  M(nx)(c^  +  k^).  Let  the  quantities 


take  their  maximum  say  for  r  =  r' ,  i  =  i' ,  j  =  j '  .  If  i'  J)r, 

the  quantity  10.6  in  the  left  member  of  10.5  may  be  replaced  by 

O'  ;  thus  we  deduce  D'  ,  <  (l+k  +  )D'  +  b*T  and 
n+1  n+1  —  o  n 


10.6  D^+1  <  (l+k0T)Dn  +  b'x  . 

If,  to  the  contrary,  i'  e  j , ,  then  by  10.4  <  Dq  :  hence 

relation  10.6  holds  in  any  event.  This  relation  also  being  true  when 
the  left  side  is  replaced  by  A^,  we  conclude  that 

10.7  D,  <  (1  +  k  x)D  +  b'T  . 

n+1  -  o  n 


Reasoning  as  at  the  end  of  section  8  proves  from  this  that 

k  nx  ,  ,  k  nt 

D  <  e  0  D  +  £-(e  0  -  1)  , 

n  -  o  k  ’ 

o 

while  D  <  max  ($>-,, A  )  .  Thus,  estimate  10.1  is  completely  justified. 
A  bound  for  rf”1 1  eru^ j  I  is  to  be  had  as  follows.  First,  apply  t)-'*' 
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to  equation  6.3  deriving  difference  equations  that  the  difference  quotients 
n""iei>u^  satisfy.  The  coefficients  in  the  difference  equations  are  bounded 
because  of  10.1  .  The  initial  values  of  the  difference  quotients  are  known, 
and  their  boundary  values  are  zero.  Then  arguing  as  in  section  8,  from  the 
difference  equations  deduce  recursive  inequalities  for 


E  =  max  tf1 1  e  uj  J  , 

n  r,i,j  r 

and,  from  the  recursive  inequalities,  relate  En  to  Eq  .  The  result  is  the 
desired  bound. 


11.  Hdlder  conditions  in  untruncated  problems.  Proof  of  Theorem  3. 

Under  hypotheses  (i— iv)  of  section  1,  any  problem  has  a  weak  solution, 

obtainable  (Section  4)  as  the  limit  of  solutions  of  truncated  problems;  if 

absolutely  continuous  with  respect  to  t,x  ,  s  =  l,...,d,  in  the  region 

ST  -  S(o)  ,  this  solution  is  strong  (Section  1).  Here,  we  shall  prove 

absolute  continuity  when  hypotheses  (i)  ,  (i).,  (iv),,  and  (v),  to  (v) ,  of 

a  l  ±4 

Section  1  are  added  to  the  others.  Under  these  added  assumptions  we  shall 
show,  more  specifically,  that  the  (weak)  solution  of  the  problem  satisfies 
a  Lipschitz  condition  with  respect  to  t  and  Holder  conditions  with  respect 
to  xg,  s  =  l,...,d,  that  sure  uniform  in  the  complement 

He  s  !ysi  >  £  >  8  =  1»***»d» 

of  an  arbitrary  e-neighborhood  of  the  planes  y  =0,  s=l,...,d.  This 
solution  accordingly  being  a  strong  solution  of  the  problem.  Theorem  3  will  be 
proved. 
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As  before,  we  normalize  c  by  the  condition  c  >  0  . 

The  Lipschitz  condition  in  t  is  immediate;  by  Section  8,  such  a 
condition  holds  uniformly  in  the  approximating  truncated  problems,  and  the 
condition  of  course  is  preserved  in  the  limit.  The  Htflder  conditions  with 
respect  to  x  remain  to  be  considered.  Since  such  conditions  also  are  pre¬ 
served  under  the  operation  of  taking  a  limit,  it  will  suffice  to  prove,  under 
the  stipulated  assumptions,  that  the  solutions  u^,T*,T(x,y,t)  of  difference 
equations  for  truncated  problems  are  subject  to  H61der  conditions  with  respect 
to  x  ,  s  =  l,...,d,  that,  for  y  e  H  and  0  <  t  <  T  ,  are  uniform  with 
respect  to  x,y,t,^,ti,'t,L.  It  will  be  enough,  in  other  words,  to  find  bounds 

for  |(X^  —  l)u^,  |/tm^)a  that  hold  for  n  t  <  T,  jr|eH  ,  and  all  £,r|,T,L. 
We  shall  do  so  vfnder  the  restriction,  for  instance,  thal  t)  =  0(£a)  . 

From  equations  6.7  we  obtain 

11.1  (1+0  x  |jsh)U“  -  Du^1  =  (l  -  txV^Kx*  -  DuJj 


+  «  JT  Bo,X=(X®  -  Du^t1  +  «  y  -  Du?*1 


sj  s  r 


ij 


5 


sj  s  '  r  ~'“ij 


+  T  x“k“  «  CX™  —  l)u“  +  T  r“  , 

1  r  ijj1  *r  IV  4  1.  -rtn  9 


|J*I<L 


r  ‘*“ij' 


ij,rm 


where 


S«,r.  -  -<?  ■  <,  ♦  -  !><,  ♦  *  E  0?  -  1W5,,.  •  « ly 

h’l<L 

Because  of  the  existence  of  the  boundary,  there  are  differences 
(X^  —  Du^1  in  right  members  of  the  above  equations  that  never  appear  on 
the  left.  These  are  of  two  types,  the  first  consisting  of  those  for  which 
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i  e|  .  ,  s  |  r;  by  the  boundary  condition,  these  vanish.  The  second  type 
s»  J 

of  difference  is  that  for  which  i  or  X  i  is  on  ^  Theorem  9.1 

r  r»  J 

estimates  this  second  kind  of  difference,  but  not  uniformly:  hence,  reason¬ 
ing  of  the  type  of  section  10  cannot  be  used  unmodified. 

For  0  <  nx  <  T,  let 

S”4m=”“i  K~lKs^  • 

where  max^  indicates  a  maximum  taken  for  such  indices  i  as  do  not 
correspond  to  a  boundary  value  (a  value  of  one  of  the  two  kinds  discussed 
above)  of  the  difference  concerned.  The  boundary  values  are  estimated  by 
Prjm  *  defined  in  section  9.  Hence,  for  the  quantity 


we  have 


Arjm  *  ■“  •  Pri«>  “h*P  *  0 


=  A.jm  rim  Jr  =  0 


|0T-l)u"|  <  <£ 


for  all  values  of  the  various  indices. 


A  uniform  estimate  for  (of  the  form  <  constant  (m^)a)  , 

according  to  a  previous  remark,  is  out  of  the  question.  For  any  e  >  0 
satisfying  condition  9*5  *  such  an  estimate  will  be  possible 9  however 9  for 


eilJrl1T£LTl  rjm 
0<|j J<  L 
(s  $  rT 
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The  process  of  estimation  will  involve,  as  auxiliary  quantities  besides  the 
latter, 


In 

rm+ 


max  n 

USI<L 

(s^r) 


I 

rm- 


3 


max  f| 

u8i<i 

(s+r) 


In 

rm 


/  Tn 

=  (lrmf  * 


Concerning  the  initial 


values  of  the  new  quantities,  it  is  pertinent  to 


state 

11.3a  E°  <  max  ({>,,  (b  +  k  M(T))/e)  ml,  , 
erm  —  J.  o  o 

as  follows  from  Theorem  9.1  and  hypothesis  (iv)^,  and 

11.3b  1°  <  ek  n£  +  C  (mOa  , 

rm  —  1  a 


a  consequence  of  Theorem  9.2. 

The  estimation  process  begins  with  the  observations  that 

|R?  |  <  c'(mOQ  , 

ij,rm 


where,  c'  =  M(T)(c  +  k  )  +  b  ,  and  that 
a  a  a 

I'^I  Oyj.  '  ■  l"dX  <  L  + 

|j'l<L  OC|jjJr£e 

(s^r) 

+  term  with  j*  =  0)|<  2k  In  +  k  e“  +  2kr|  M(T)  . 
r  *—  rm  o  erm 


z 

l^h>e 
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From  these  and  11.2  ,  equation  11.1  easily  gives  us  the  preliminary  in¬ 
equality 


(1+6  T  |j  h)  a1*1  <  An  |j  h  +  2kxl"  +  2xr)k  M(T) 

'  '"s'  -rim  —  rim  '"s'  1  r.lm  rm  1 


's''1'  rjm  —  rjm  ~  ^  '•V  rjm 


+  k  xEn  +  c'T(mOQ  • 
o  erm 

If,  when  jy  ^  0.  S  be  replaced  by  P^m  ,  the  inequality  remains 

good.  Hence,  A  can  be  replaced  by  A^^  .  Making  the  replacement  and 
a  consequent  cancellation,  we  arrive  at 

11.4  A “t*  <  A".  +  x(2k  l“  +  k  Ep  +  c'  (m^)a  +2r,  k  M(T) )  . 

rjm  —  rjm  rm  o  erm 


Summing  appropriately  over  j  ,  we  now  easily  deduce 


Ci  <  +  ‘ C, +  Ek„  +  e<=’^)a+  2i-'k  M<T» 


and  a  similar  inequality  concerning  I*  ;  from  both,  we  have 


11.5  I**1  <  (l+2kre)l”  +  xe(k  e”  +  C  0*)°  +  2qEk  M(T) )  . 

rm  —  rm  o  Erm  ’ 


From  11.4  we  also  have 


11.6  E^  <  (l+H)Enm  +  T(2k  l“  +  c'(mOa  +  2tje  k  M(T)) 
erm  —  o  erm  rm  '  '  ' ' 


Hence,  for 


F  =  Fn  =  max  (ln  ,  En  ) 
n  rm  rm'  srra 


we  can  deduce  the  recursive  inequalities 
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Fn+1  <  (1  +  ai:)Fn  +  T[c'(mOa  +  2tiek  M(T)J  , 

where  a  =  (kQ  +  2k) max  (l,e)  .  By  an  inductive  process  like  that  of 

section  8, 

anT  , 

Fn  <  eanTFo  +  (® - J-i)  [c‘(mOa  +2r,ek  M(T)]  , 

while  Fq  can  be  estimated  from  11.3  .  The  result  is  a  bound  for  Fn/(mOa, 
and  thus  for  A*V/(mS)°  ,  depending  only  on  the  constants  of  the  problem. 
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NOTES 

[p.37]  3/  See  the  discussion  and  references  in  Section  6,  Douglis, 
On  discontinuous  solution  of  quasi-linear,  first  order 
partial  differential  equations,  NAVORD  Report  6775,  U.  S. 
Naval  Ordnance  Lab.,  White  Oak,  Md.,  22  Jan.,  1960. 
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